finish
/clear
/filn, E323  
/title, FEA of connecting zone of support of evaperator
! ********* 趨*********
Dci=5000                              !Ͳھ
tc=40                                 !Ͳں
tcon=40                               !׶ηͷں
alfa=30                               !׶ηͷ׶
lc=2500                               !Ͳ峤

*afun,deg                             !Ƕȵλ趨
Lcon=Dci/2/tan(alfa)                  !׶ηͷȫ׶
lrat=0.5                              !ģ׶ηͷ߷򳤶ϵ
Lcon=lcon*lrat                        !ģ׶ηͷ߷򳤶

ns=4                                  !֧
h0=1500                               !ֱ֧װλ
Dis=Dci+tc                            !֧ˮƽװλ
ts=20                                 !֧塢
ws=500                                !
hs=1300                               !߶
bb=800                                !װ峤
w1=250                                !λ
pi=0.1                                !ѹֵ
Gra=9.81                              !ٶȣm/s2
denstell=7850                         !ֲܶȣkg/m3
denwater=1000                         !ˮܶȣkg/m3
m1=58.1e3                             !أkg
hy1max=13                             !׶ηͷƽҺ߶
hy0max=15.0                           !Һ߶
py1=hy1max*denwater*Gra/1e6           !׶ηͷƽҺѹ
py0=hy0max*denwater*Gra/1e6           !Һѹ
scon=(Dci/2*lrat*2+tcon*cos(alfa))*tcon	!ģ׶ηͷ/?
pcon=(pi+py0)*(Dci/2*lrat)**2/scon/cos(alfa)            !ģ׶ηͷƽغ
pcm1=m1*gra/(3.1415927*((Dci+2*tcon)**2-Dci**2))        !Ͳغ
pc=pi*Dci**2/((Dci+2*tcon)**2-Dci**2)-pcm1              !Ͳƽغ
!****************ǰ***************************
/prep7	
et,1,45                               !嵥Ԫ
mp,ex,1,2e5                           !ϵĵģ
mp,nuxy,1,0.3                         !ϵĲɱ
!****************ģ***************************
k,1,0,0,0                             !תĹؼ
k,2,0,Lc,0                            !תĹؼ
BLC4,Dci/2,0,tcon,-lcon               !׶ηͷĸ
local,11,1,Dci/2                      !ֲϵ
AGEN,2,1,,,,-alfa,,,,1                !תĸõת׶ηͷ
BLC4,Dci/2,0,tcon,lc                  !Ͳ
NUMMRG,ALL, , , ,LOW                  !ϲͬ
ksel,s,loc,y,-tcon,0                  !ѡͲͷؼ
*get,k1,kp,,num,min                   !ȡؼ
ksel,u,,,k1
*get,k2,kp,,num,min
*get,k3,kp,,num,max
ksel,a,,,k1
a,k1,k2,k3                            !ɺ
alls                                  !ȫѡ
vrotat,all,,,,,,1,2,360/2/ns          !Ͳ塢ͷ1/8ת
wpoff,0,0,-bb/2                       !ƶ
BLC4,,-h0,Dis/2,hs,ts                 !ɽĸ
BLC4,Dis/2-w1,-h0,-ts,hs,bb/2         !ɸĸ
BLC4,,-h0,Dis/2,ts,bb/2               !ɵװĸ
vsel,s,,,4,6                          !ѡ塢弰װĸ
vsba,all,5                            !׶ηͷи塢弰װ
vdele,7,,,1                           !ɾĸಿ
vdele,10,,,1                          !ɾĸಿ
vdele,11,,,1                          !ɾװĸಿ
wpoff,Dis/2-ws                        !ƶ
wprot,,,90                            !ת
vsbw,all                              !ٴиװ
vdele,4,,,1                           !ٴɾಿ
vdele,6,,,1                           !ٴɾװಿ
vptn,all                              !塢弰װ廥
alls
!*****************************************
KWPAVE,23                             !ƽؼ23
wprot,,,90                            !תƽ
vsbw,all                              !ùƽʷ
KWPAVE,19                             !ƽؼ19
vsbw,all                              !ùƽʷ 
wprot,,90,                            !תƽ 
vsbw,all                              !ùƽʷ 
vptn,all                              !廥
alls                                  !ȫѡ
lglue,all                             !ճϸ
lesize,36,,,3,,,,,                    !趨߷ֶ
lesize,55,,,3,,,,,                    !趨߷ֶ
lesize,45,,,3,,,,,                    !趨߷ֶ
esize,40                              !趨ԪС
vsweep,all                            !ӳ仮
nummrg,all                            !ϲ
numcmp,all                            !ѹ
finish
!*******************************************
/solu
wpcsys,-1,0                           !תԭʼϵ
asel,s,loc,z,0                        !ѡz=0
da,all,SYMM                           !ʩӶԳԼ
local,12,1,,,,,-90                    !ֲϵ
asel,s,loc,y,360/ns/2                 !ѡ1/8λý
da,all,SYMM                           !ʩӶԳԼ
csys,0                                !ֱϵ
asel,s,loc,y,-h0                      !ѡ֧
asel,r,loc,x,Dis/2-ws,Dis/2
da,all,all                            !Լλ
asel,s,loc,y,Lc                       !ѡͲ
sfa,all,1,pres,-pc                    !ʩƽغ
asel,s,loc,y,-Lcon*cos(alfa)-tcon,-Lcon*cos(alfa)  !ѡ׶ηͷ
sfa,all,1,pres,-pcon                  !ʩƽغ

asel,s,area,,83,87,4                  !ѡͲڱ
asel,a,area,,42,
sfa,all,1,pres,pi                     !Ͳڱʩѹ
asel,s,area,,44,46,2,1                !ѡͷڱ
asel,a,area,,93,95,1
asel,a,area,,91
sfa,all,1,pres,pi+py1                 !׶ηͷڱʩѹҺѹ
alls
solve                                 !
!*******************************************
/post1
PLNSOL,S,INT,0,1                      !ʾӦͼ
fini
